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We present a search for the rare B-meson decay B° —* afp^ 1 with af — > ir + n ir . We use (110 ± 
1.2) x 10 6 T(4S) -» BB decays collected with the BABAR detector at the PEP-II asymmetric-energy 
B Factory at SLAC. We obtain an upper limit of 30xl0~ 6 (90% C.L.) for the branching fraction 
product B{B° -> af p^)B(af — > 7T+7T ty ± ), where we assume that the af decays exclusively to 
p°n ± . 



PACS numbers: 13.25.Hw, 12.15.Hh, 11.30.Er 

In the Standard Model, CP-violating effects in the B- 
meson system arise from a single phase in the Cabibbo- 
Kobayashi-Maskawa (CKM) quark- mixing matrix [1]. 
The decay B° — > afp T proceeds via a b — > uud tran- 
sition [2], and interference between direct decay and 
decay after B°B° mixing results in a time-dependent 
decay-rate asymmetry that is sensitive to the angle a = 
8,rg[-VtdV* b /V ud V* b } [3] in the unitarity triangle of the 
CKM matrix. An additional motivation for studying 
B a — > af p T is that this is a significant background to 
B — > p/9 decays, e.g. [4-8], which currently provide the 
most accurate measurement of a. The ARGUS exper- 
iment previously searched for the decay B° — > afp^, 
which resulted in an upper limit of B(B° — > p T ) < 
3.4 x 10" 3 (90% C.L.) [9]. This paper presents the result 
of a search for B° — > af with af — > Tr" 1 "^ - ?^, where 
we assume that the af decays exclusively to / 9°7r ± . A 
theoretical prediction of the branching fraction B(B° — ► 
afp*)B{af — > (37r) ± ) has been made by Bauer, Stech 
and Wirbel [10] within the framework of the factoriza- 
tion model. They predict a value of 43xl0~ 6 , assuming 
\V ub /V cb \ = 0.08. 

The data used in this analysis were collected with the 
BABAR detector at the PEP-II asymmetric-energy B Fac- 
tory at SLAC during the years 2003-2004. This repre- 
sents a total integrated luminosity of 100 fb" 1 taken at 
the T(4S) resonance (on-peak), corresponding to a sam- 
ple of 110 ± 1.2 million BB pairs. An additional 21.6 
fb -1 of data, collected at approximately 40 McV below 
the T(4S) resonance (off-peak), were used to study back- 
ground from e + e~ —tqq{q = u, d, s, c) continuum events. 

The BABAR detector is described in detail else- 
where [11]. Surrounding the interaction point is a sil- 
icon vertex tracker (SVT) with 5 double-sided layers 
which measures the impact parameters of charged par- 
ticle tracks in both the plane transverse to, and along 
the beam direction. A 40 layer drift chamber (DCH) sur- 
rounds the SVT and provides measurements of the trans- 
verse momenta for charged particles. Both the SVT and 
the DCH operate in the magnetic field of a 1.5 T solenoid. 
Charged hadron identification is achieved through mea- 
surements of particle energy-loss in the tracking system 
and the Cherenkov angle obtained from a detector of 
internally reflected Cherenkov light. A CsI(Tl) electro- 
magnetic calorimeter (EMC) provides photon detection, 
electron identification, and 7r° reconstruction. Finally, 
the instrumented flux return of the magnet allows dis- 
crimination of muons from pions. 



We reconstruct B° — ► af p~ candidates from combi- 
nations of af — ► 7r + 7r~7r + and p~ — ► 7r°7r~ candidates. 
The ai(1260) — > 37r decay proceeds mainly through the 
intermediate states {-ktt) p -k and (iti:)^ [12]. We do not 
distinguish between the dominant P-wave (7T7r) p and S- 
wave (7^)0- in the channel 7r + 7r~. The Monte Carlo (MC) 
signal events are simulated as B° decays to af (1260)p~ 
with af -> p°ir+ using the GEANT4-based [13] BABAR 
MC simulation. Possible contributions from B° decays 
to a^(1320)p~ and 7r + (1300)p~ are investigated. 

We only consider events that have a minimum of one 
7T° and four charged tracks, where the charged tracks are 
required to be inconsistent with lepton, proton and kaon 
hypotheses. 

We form 7r° — > 77 candidates from pairs of photon 
candidates that have been identified as localized energy 
deposits in the EMC that have the lateral energy distri- 
bution expected for a photon. Each photon is required 
to have an energy E 1 > 50 MeV, and the tt° is required 
to have an invariant mass of 0.10 < m 77 < 0.16 GeV/c 2 . 

The p~ mesons are formed from one track that is con- 
sistent with a 7T~ and the aforementioned tt° candidate. 
The candidate p~ is required to have an invariant mass 
of 0.5 < m p - < 1.1 GeV/c 2 . We also constrain the cosine 
of the angle between the tt° momentum and the direction 
opposite to the B° in the p~ rest frame (cos# p -) to be 
between —0.9 and 0.98. This removes backgrounds which 
peak at the extremes of the distribution where the signal 
reconstruction efficiency also falls off. 

We form the af candidate from combinations of three 
charged pions. We first form a p° — > 7r + 7r~ candidate 
from two oppositely charged tracks. This combination 
is required to have an invariant mass of 0.4 < m p o < 
1.1 GeV/c 2 . The af candidate is then formed by adding 
another charged track to the p°, and requiring that the 
mass of the af satisfies 0.6 < m a + < 1.5GeV/c 2 . The 
vertex of the B-candidate is constrained to originate from 
the beam spot. In order to reduce background from con- 
tinuum events we require that | cos(#t)| < 0.7, where 9t 
is the angle between the B thrust axis and that of the 
rest of the event (ROE). 

We use two kinematic variables, wes and AE, in order 
to isolate any signal. We define the beam-energy substi- 
tuted mass toes = \A\/*/2) 2 ~ (p^) 2 , where s/s is the 
e + e~ center-of-mass (CM) energy. The second kinematic 
variable, AE, is the difference between the B-candidate 
energy and the beam energy in the CM frame. We re- 
quire toes > 5.25GeV/c 2 and -0.15 < AB < 0.1 GeV. 
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Additional separation between signal and continuum 
is obtained by combining several kinematic and topolog- 
ical variables into a Fisher discriminant (J 7 ) [14]. The 
variables L , L 2 , and | cos 8tr\, and the output of a mul- 
tivariate tagging algorithm [15] are used as inputs to T . 
L and L 2 are defined as 

L =5>il. £ 2 =]>>*|cos(^) 2 , (1) 

ROE ROE 

where the sum is over the ROE, p* is the particle mo- 
mentum in the CM frame. 6i is the angle of the particle 
direction relative to the thrust axis of the i3-candidate, 
and cos 9tr is the cosine of the angle between the B 
thrust axis and the beam axis. The multivariate tagging 
algorithm identifies the flavor of the other B in the event 
to be either a B° or B°. The output of this algorithm is 
ranked into categories of different signal purity. 

We expect the polarization of the a\p~ final state to be 
predominantly longitudinal, as was found in the similar 
decay B — > pp [4-8]. We have used both longitudinal 
and transverse polarized signal MC simulated data in 
this analysis. After applying the selection cuts above, we 
have 2.8 (2.3) longitudinal (transverse) polarized signal 
MC simulated data candidates per event. 

We define as self-cross-feed (SCF) the set of candi- 
dates that were incorrectly reconstructed from particles 
in events that contain a true signal candidate. We select 
one B candidate per event in which the mass of the re- 
constructed p° is closest to that of the true p° mass [16]. 
Choosing the candidate using the p° mass reduces the 
SCF fraction by 18% relative to a random selection. To 
avoid potentially biasing our final result, we do not use 
information from the p° meson in the remainder of the 
analysis. After all selection cuts have been applied, the 
longitudinal and transverse signal SCF fractions are 0.58 
and 0.42, respectively. The selection efficiency of longi- 
tudinal (transverse) signal is 9.44% (10.15%). 

Besides the continuum background we also have back- 
ground from B decays. We divide the B-background 
into the following four categories according to £>-mcson 
charge and the charm content of the final states: (i) B° — ► 
charm, (ii) B° — > charmless, (iii) B^ — > charm and (iv) 
B ± — ► charmless. From large samples of inclusive MC 
simulated data we expect 2394, 424, 3281 and 215 events 
of these background types, respectively. In addition, a 
number of exclusive _B-background modes that have a 
similar final state to the signal were studied. This in- 
cludes those that have an intermediate a\ meson in the 
decay. None of these modes were seen to have a signifi- 
cant efficiency after the selection cuts had been applied. 

We perform an extended unbinned maximum likeli- 
hood fit to the data. The likelihood model has the 
following types: (i)-(iv) the four aforementioned in- 
clusive _B-background categories, (v) true signal, (vi) 
SCF signal and the (vii) e + e~ — > qq (q = u, d, s, c) 



continuum background. The probability density func- 
tion (PDF) for each event i has the form Pj )C = 
-Pi. c ( TO ES7 Ai^-7 r ! m a +! m / !)-! cos #p-)- From these indi- 
vidual PDFs the total likelihood 

n 

L = e-"'H^«,c, (2) 

i—l c 

is constructed. The parameters n and n' are the numbers 
of selected on-peak events, and the sum of the yields N c , 
where c is one of the seven types in the likelihood model. 
Correlations between the m p - and cos 9 p - variables are 
taken into account for the real (T) and fake combinato- 
rial (F) p~ candidates. All other correlations between 
fit variables are seen to be small. However, the effect of 
ignoring them results in a bias on the fitted signal yield 
which is discussed below. 

Each of the signal distributions has a signal yield and 
a polarization fraction that are denoted by N S i g (N' sig ) 
and /l (/l), for the true (SCF) signal, such that the sum 
of the true and SCF signal is described by 

N sig [f L P<° n9 > true + (1 - /L )pf ""<<"-] + 
N' 8tg [f L P!° n9 ' SCF + (1 - f L )P? an ' SCF }. (3) 

The continuum yield, N S i g , N' sigl and the parameters of 
the continuum rriES and AE PDFs are allowed to vary in 
the fit. Under the assumption that no significant signal 
is observed, the value of Jl is fixed to 1.0 in the fit. The 
value of f' L is also fixed to 1.0 in the fit since it is highly 
correlated with N S i g , N' sig and Jl ■ Only the fitted value 
of N sig is used to derive the final result. We also fix the 
B-background yields to the aforementioned values. 

The PDFs used for each component are given in Ta- 
ble I. The signal and ^-backgrounds are parameterized 
using MC. We use a non-parametric smoothing algo- 
rithm [17] when defining some of the background PDFs 
(as indicated by the abbreviation NP). We account for 
the difference between F and T p~ — ► tt~tt° distribu- 
tions in the background using the product of ID PDFs, 
denoted in Table I as 'BG m-hel', such that 

Pi (m p -, cos dp-) = (1 - a T )PF,i{m p -)PF,i(cos6 p -) + 
a T Pr,i (m p - )P T ,i (cos 9 p - ) , (4) 

where ax is the fraction of T events. The continuum 
shape for cos6> p - {m p -) is derived from off-peak (on- 
peak) data. The true p~ resonance Breit Wigner shape 
uses m p - = 0.77 GcV/c 2 , and V = 0.150 GeV [16]. The 
parameterizations used for this PDF are summarized in 
Table II. 

The results from the fit are N slg = 90±38(stat), N' sig = 
42 ±98 and a continuum yield of 25798 ± 182 events. The 
bias on the fitted signal yield is evaluated by perform- 
ing ensembles of mock experiments using signal MC em- 
bedded into MC samples of background generated from 
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TABLE I: The types of PDFs used to model the different variables for each component in the likelihood fit, where the PDFs 
underlined have their parameters varying in the nominal fit. The abbreviations are: G = Gaussian, G2 = Double Gaussian, 
G3 = Triple Gaussian, CB = Crystal Ball (a Gaussian with a low side exponential tail) [18], ARGUS = ARGUS function 
— x 2 exp [— £(1 — x 2 )] , with x = 2m-Es/yfs ancl parameter £ [19], which is allowed to vary in the fit, Pn = Polynomial 
of order n, BW = Breit-Wigner, helicity = cos 2 - or sin 2 6 p - depending on partial wave which is modified by a quadratic 
acceptance function, BG m-hel = Background cos 6 p - and PDF of Eq. 4, off-peak = PDF taken from off-peak data, and 
ID = smoothed ID histogram. 



Component mES AB T m a + cos 9 p - m p ~ 



signal (long/trans. /true/SCF) CB CB+G G2 G3 helicity BW+P4 

qq ARGUS PI G2 (off-peak) ID (off-peak) BG m-hel BG m-hel 

B° (B+H charm (charmless) NP NP NP NP BG m-hel BG m-hel 



TABLE II: The types of PDFs used to model the different 
background cos - and m p - PDF shapes. The abbreviations 
Pn and BW are denned in the caption of Table I. 



Component 



T cos e„ 



F cos6> - F m„ 



qq 




P2 


BW + PI 


P5 


P4 


B ± ^ 


charmless 


P4 


BW 


P5 


P3 


B° 


charmless 


P4 


BW 


P5 


P3 


B ± ^ 


charm 


P2 


BW 


P5 


P3 


B° -> 


charm 


P2 


BW 


P5 


P3 



the PDF. The bias is found to be ±22 events (24%), re- 
sulting in a corrected signal yield of 68 ± 38(stat). In 
Fig. 1 we compare the true signal and continuum PDF 
shapes (solid curves) to the data (points) using the event- 
weighting technique described in Rcf . [20] . The distribu- 
tions shown in Fig. 1 are not corrected for fit bias, and 
the uncertainty on each of the data points is statistical. 
No change in signal yield is seen when a^(1320)/9~ and 
7r + (1300)p~ components are included in the fit. 

Table III summarises the systematic uncertainties on 
the signal yield. Each entry in the table indicates one sys- 
tematic effect, and the contributions are added in quadra- 
ture to give the total presented. The uncertainty due to 
PDF parameterisation is evaluated by variation of both 
the signal and the background PDF parameters within 
their uncertainties about their nominal values. The un- 
certainty from the continuum tues and AE PDFs which 
are allowed to vary in the fit, are only included in the 
quoted statistical uncertainty. We assign a systematic 
uncertainty due to fit bias, evaluated as half of the fit bias 
correction on the signal yield. To validate the expected 
-B-background yields, and to assign a systematic uncer- 
tainty we perform a number of cross-checks in which we 
allow the background yields to vary in turn when fitting 
the data. We use a control sample of B — > Dp events to 
determine the systematic uncertainty in the fraction of 



SCF signal events. The effect of exclusive B meson de- 
cays to final states including ai-mesons were evaluated 
using ensembles of mock experiments. In particular, the 
systematic uncertainty on the signal yield from neglect- 
ing B — > aiai modes in the fit is 6 events. We assign 
a systematic uncertainty from using a relativistic Breit- 
Wigner with a Blatt-Weisskopf form factor with a range 
parameter of 3.0 GeV -1 for the a+ meson line shape. 
In the fit we assume that the a+ meson width, T a +, is 
400 MeV. We evaluate a systematic uncertainty due to 
this assumption by varying T a + over the experimentally 
allowed range: 250 - 600 MeV [16]. The difference in 
the distribution of T between data and MC is evaluated 
with a large sample of B — > D* p decays. The systematic 



TABLE III: The systematic uncertainties on N S i g (events). 



Source 



Uncertainty on N S[ 



PDF parameterisation 
Fit bias 

B-background yields 
SCF fraction 

Neglecting B — > cnai modes in fit 
af line shape 
af width 

Fisher data/MC comparison 



+27 
-30 
±11 
+29 
-42 
±7 
±6 

±10 

±9 
±6 



Total 



+45 
-56 



uncertainties that contribute to the branching fraction 
only through the efficiency come from charged particle 
identification (6.0%), tt° meson reconstruction (3.0%), 
tracking efficiency (3.2%), and the number of B meson 
pairs (1.1%). The systematic error contribution from MC 
statistics is negligible. 

When the fit bias correction of —22 events is ap- 
plied to the signal yield, and one accounts for sys- 
tematic uncertainties, the significance of the result is 
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FIG. 1: The true signal (top) and continuum (bottom) distributions for (left to right) tues, Ai5, m a + using the weighting 
technique described in Ref. [20]. The points represent the weighted data, and solid curves represent the corresponding PDFs. 



0.95 standard deviations. Figure 2 shows the distribu- 
tion of — ln(L/L max ) for the fit, with and without these 
systematic errors. L max is the value of the likelihood 
corresponding to the nominal fit result. The branch- 
ing fraction value for the fit-bias-corrected signal yield 
of 68 ± 38(stat)+5g(syst) is B{B a -» a+p-)B(a+ -> 

7r+7r-7r+) = (15.7±8.7(stat)+^-g(syst))xl0- 6 . This as- 
sumes that Jl — 1.0 and that the branching fraction 
of al — > n + n~n + = 0.5. As the signal yield obtained 
is not significant, we calculate the upper limit xtjl, by 




Signal Yield 



FIG. 2: The — ln(L/L max ) distribution from the fit to data 
with /_l = 1.0 (fixed). This distribution has been corrected for 
fit bias. The solid curve is for statistical errors only, and the 
dashed curve includes systematic errors. 



integrating the likelihood function (including systematic 
uncertainties) from to xtjl, for different physically al- 
lowed values of such that the C.L. of the upper limit 
is 90%. As the signal efficiency is a function of /j,, we 
report the most conservative upper limit obtained, which 
corresponds to ft, = 1.0. On doing this, an upper limit 
of 30xl0~ 6 (90% C.L.) is obtained. 

We have performed a search for the decay B° — ► afp T 
in a data sample of 100 fb^" 1 . After correcting for fit bias 
and accounting for systematic uncertainties, the signal 
yield is 68 ± 38(stat)^g(syst) events, with a significance 
of 0.95cr. As there is no significant evidence for a sig- 
nal, we place an upper limit of 30xl0 -6 (90% C.L.) on 
B(B° — > a~l p~)B{a~l — > 7r + 7r _ 7r + ), where we assume that 
the decays exclusively to p Q n + . Assuming B(a^ — > 
7r + 7r~7r + ) is equal to B(a^ — > 7r + 7r°7r°), we obtain 
B(B° -> a+p-)B{a+ -> (3tt)+) < 61xl0- 6 (90% C.L.). 
This upper limit corresponds to a significant improve- 
ment over the previous bound and is compatible with 
theoretical expectations [10]. This result is a significant 
improvement in constraining an important B background 
contribution in B — > pp decays. 
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